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Abstract 



Using the density matrix renormalization group, we study domain wall 

structures in the t-J model at a hole doping of x = |. We find that the 

domain walls are composed of d x 2_ y 2 pairs and that the regions between the 

domain walls have antiferromagnetic correlations that are ir phase shifted 

across a domain wall. At x = g, the hole filling corresponds to one hole per 

two domain wall unit cells. When the pairs in a domain wall are pinned by an 

external field, the d x 2_ y 2 pairing response is suppressed, but when the pinning 

is weakened, d x 2_ y 2 pair-field correlations can develop. 
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In the low temperature tetragonal (LTT) phase of Lai.6- a; Ndo.4Sr a ;Cu04, the tilt pattern 
of the CuOg octahedra form lines of displaced oxygens parallel to the Cu-0 bond directions. 
These lines are rotated by 90° between adjacent layers. At a filling of x = |, superconductiv- 
ity is suppressed and neutron scattering studies |1],0 reveal a striped domain wall ordering 
of holes and spins which is believed to be commensurately locked by the tilt distortion of 
the lattice. One model for this striped order |l|J^] is illustrated in Fig. 1(a). Here the charge 
domain walls are shown running vertically and centered along the Cu-O-Cu legs, although 
the phase information required to determine whether the domains should be leg centered 
or bond centered (centered between two legs) is not known. As shown, the domains are 
separated by four Cu-O-Cu spacings and for x — | contain one hole per two 4x1 domain 
wall unit cells. This latter feature is at odds with one-electron Hartree-Fock calculations 
|3|-[8|] which predict a domain wall filling of one hole per domain wall unit cell. The spins 
in the regions between the walls are antiferromagnetically correlated with a ir phase shift 
across a domain wall. When x ^ |, superconductivity is found to coexist with a weakened 
domain wall ordering, suggesting a close connection between the two. 

Here we present results of numerical density matrix renormalization group calculations 
for at- J model with a hole doping x = |. We find evidence for domain walls with tt phase- 
shifted antiferromagnetic regions separating the walls, and with a filling of one hole per two 
4x1 domain wall unit cells. In contrast to Fig. 1(a), however, the domain walls are bond 
centered, and are made up of d x 2_ y 2 pairs of holes. We find that just as for the two-leg 
and four-leg ladder problems, |9"1,|T0|1 there is a competition between a charge density wave 
(CDW) phase and a superconducting pair phase. The Coulomb interaction, arising from the 
90° rotated domain walls on adjacent planes in the LTT phase PJUfl, acts to pin the pairs 



within a domain wall for x — |. We show that when this pinning is weakened, d x 2_ y 2 pair 
field correlations develop. 

The t-J Hamiltonian we will study is given by 

H = -J2 tp c (cl, s Cj,s + c),sd,s) P G + J2 J 4 - J n i n j) ■ (!) 

(U)b (ij) V 4 J 



Here (ij) are near-neighbor sites, s is a spin index, Si = c\ s <T S:S >c ijS > and = c\^c^ + cj^c^, 
with c\ s (q s ) an operator which creates (destroys) an electron at site % with spin s. The 
near-neighbor hopping interaction is t, the near-neighbor exchange interaction is J, and the 
Gutzwiller projection operator Pq excludes configurations with doubly occupied sites. We 
refer to the Cu-Cu lattice spacing as a and measure energies in units of t. 



Our calculations were carried out using a density matrix renormalization group [12 



(DMRG) algorithm in which transformation matrices were stored and used to construct the 



initial state for each superblock diagonalization ||13|| . We present results here for 8 x 8 and 
4 x 18 doped systems. For the more difficult 8x8 system, typically 10 equilibration sweeps 
were made through the lattice, keeping of order 10 3 states per block on the final sweep. The 
transformation matrices were used to calculate the ground state expectation values of the 
desired operators at the end of the calculation. DMRG is extremely accurate for narrow sys- 
tems, but its accuracy decreases rapidly as the width increases, and consequently, our 8x8 
systems may represent the most ambitious DMRG calculations to date. Truncation errors 
were approximately 0.0003, which we consider to be at the limits of acceptability. In addi- 
tion, the calculations sometimes got stuck in metastable hole configurations, in cases where 
two or more states had a very close separation in energy but corresponded to completely 
different hole configurations. In these cases it was necessary to perform several calculations, 
with the holes initially localized in different configurations, using the total energy to choose 
between different runs. (In contrast, the spin configuration equilibrated more readily.) Nev- 
ertheless, the domain wall configurations shown here were robust as long as the initial hole 
configurations were reasonably close to the final configurations displayed. In particular, 
three of the most robust features emerging from a variety of runs were that the domain 
walls are bond centered, composed of d x 2_ y 2 pairs, with it phase shifted antiferromagnetic 
regions separating the walls. 

Figure 1(b) shows the charge and spin density in the ground state of the 8x8 system 
with J/t = 0.35 and 8 holes, corresponding to a filling x — |. Periodic boundary conditions 
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were used in the y-direction, and open boundary conditions in the x-direction. Along the 
left and right edges of the system a small staggered magnetic field of 0.05t was applied. 
The boundary conditions and the edge staggered field serve to orient and pin the domain 
walls in the configuration shown. In an LTT phase, the increased Cu-0 bond length in one 
direction would act to orient the domain walls through an anisotropic hopping and exchange, 
ty = jtx and J y = 7 2 </r, with 7 « 1, but for simplicity we have used 7 = 1 throughout. 
The staggered edge field further acts to pick a direction for the spin order, which allows 
direct measurement of the spin configurations and reduces truncation errors in the DMRG 
calculation. The charge density in the x-direction shows a strong modulation with period 
4a, and the filling of the wall is consistent with that shown in Fig. 1(a). In addition, a 
small charge density modulation, also with period of 4a, is present in the y direction ||. 
This modulation is pinned in the y direction despite periodic boundary conditions by the 
truncation errors in the calculation. We cannot determine whether an exact calculation of 
the two-dimensional t-J lattice would have small pinned charge density modulations along 
the walls, but it seems clear that the system is near a CDW instability, with period 4a, 
along the walls. The spin response (Sf) is also shown in Fig. 1(b) and corresponds to 
antiferromagnetic correlations which have a 7r phase shift across the bond domain wall. 
This 7r phase shift was purely the result of the simulation, not enhanced in any way by 
boundary or initial conditions. 

In Fig. 1(c), the black circles show the most probable configuration of the eight holes 
obtained by maximizing the ground state expectation value of 

p(< 1 ,4,...,40 = ILp(4), ( 2 ) 

i 

with p(£) = (1 — n^)(l — n^) the hole projection operator for the £ th lattice site. The 
thickness of the lines connecting various sites denotes the strength of the exchange field 
(Si ■ Sj) when the holes occupy the most probable configuration, consisting of four pairs. 
Note that the most probable configuration of a pair is a diagonal configuration, with a 
strong exchange bond running diagonally between the spins adjacent to the holes. This 
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configuration is the most likely configuration of a pair in a variety of t-J clusters [|14|,|10 



From a strong coupling point of view, the pairing can be viewed as arising from a com- 
promise in which two holes locally arrange themselves so as to minimize the disturbance of 
the background exchange energy while at the same time lowering their kinetic energy. In 
the unphysical regime of J > t, the holes would tend to be near neighbors to reduce the 
number of broken exchange bonds. In the physical region J < t, the kinetic energy plays 
an increasingly important role so that, as shown in Fig. 1(c), the most probable pair config- 
uration has the holes sitting on diagonal sites with strong singlet correlations on the other 
diagonal. In this case, four of the eight one-electron hops leads to a configuration with a 



near-neighbor exchange bond which stabilizes the pair. As previously discussed, fT5|J14]| the 
pair structure is such that it has an overlap with the undoped antiferromagnetic background 
through a hole pair field operator which has d x 2_ y 2 symmetry. Here we will see the d x 2_ y 2 
structure of the pairs from the pair field response discussed below. 

In the LTT phase, the CuOg tilt structure causes the domain walls to be perpendicular 
in adjacent planes [|TJ. This gives rise to an electrostatic potential, with a period of 4a, along 
the domain wall [II|. For x = 1/8, this corresponds to the CDW instability of the d x 2_ y 2 
domain wall and can lead to a pinning of the pairs along the wall. Here we have modeled 
this effect by adding a spatially varying site potential AV = O.lt on the 8x8 lattice shown 
in Fig. 2. The sites with the extra potential AV are shown by the shaded rectangles, which 
also indicate the domain walls in the adjacent planes. The potential pins the pairs, forming 
a CDW lattice of pairs. 

The charge and spin structure factors 

s c(q) = ^E e<q/ W + n ii) ( 3 ) 

and 

SM) = ^Ee^(n tl -n ri ) (4) 

for this lattice are also shown in Fig. 2. The intensity in S c (q) at (■£-, 0) is expected for ver- 
tical charged domain walls. The intensity at (j-, j-) and (0, reflect the pair correlations 
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along the domain walls. With the 90° rotations of the domain walls from layer to layer, it is 
the ( j-, y) peak that would be important to observe. The spin structure factor at (§^, f ) 
reflects structure of the 7r-phase shifted antiferromagnetic regions. The amplitude of the 
third harmonic intensities is much weaker. 

At dopings away from x = 1/8, the Lai. 6 _ x Ndo.4Sr x Cu04 system becomes superconduct- 
ing, and the size of the tilt modulation as well as the intensity of the magnetic Bragg peaks 
decreases [|J. This suggests that when the pinning is weakened, either through a reduced 
pinning potential or through a mismatch in the periods of the pinning potential and the 
CDW instability, stripe order and superconducting pairing can coexist. As a test of this, we 
have used a 4 x 18 cluster to model a single, longer domain wall, with a variable pinning 
potential AV which acts on the sites in the shaded regions shown in Fig. 3(a). This cluster 
has open boundary conditions on all sides with a staggered magnetic field of magnitude 0.05£ 
applied along the top and bottom edges. The magnetic field has a 7r phase shift between 
the edges in order to mimic the single domain wall structure shown in Fig. 2(a). In order 
to measure the tendency for superconductivity, a weak pair field proximity effect term 

H x = dJ2(A\. +§ + A^) (5) 

i 

was added to the Hamiltonian. Here 

M,i+y — C i,1 C i+y,l + C i+yA C i,l (^) 

adds a singlet electron pair to sites i and i + y. Note that the interaction H\ couples equally 
to d x 2_ y 2-Y\ke and extended s-wave-like pairs. That is, it does not distinguish whether the 
internal structure of the pair field has a change in sign for the singlet components in the 
y-direction relative to those in the x-direction. In order to include this term, rather than 
use the number of electrons N as a good quantum number to break up the Hilbert space, 
N modulo 2 was used. (Total S z was conserved in the usual fashion.) We then measured 
the strength of the pair fields in the ground state in both the x and y directions, (A^+a) 
and (Ai t i + y) for all sites i. The charge and spin structure of the 4 x 18 cluster with 8 holes 
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is shown in Fig. 3(a) for A V = 0.05t. In Fig. 3(b), the pair field strength is shown by the 
width of the lines, and the sign of the field is indicated by the type of line, dashed or solid. A 
relative sign difference between the x and y directions indicates d x 2_ y 2 pairing. In Fig. 3(b) 
we see a significant d x 2_ y 2 pair response coexisting with a modest charge density wave. In 
Fig. 3(c), we show the pair field (A m y) averaged over all the ?/-bonds in the middle four 
rungs as a function of AV^. The suppression of pairing by the CDW is evident. We expect 
that a larger, two dimensional array of domain walls would show a more enhanced response 
versus AV. 

These calculations show that holes doped into at- J lattice can form domain walls of pairs. 
For a filling x = |, these walls have an average filling of one hole per two domain wall unit 
cells, and there is a tendency for the pairs to form a pinned CDW structure. If the pinning 
is weakened, the pairs fluctuate, developing phase coherence, and d x 2_ y 2 superconducting 
correlations appear. The idea that the domain walls are made up of pairs provides a natural 
explanation of the special filling x — | and the intimate relationship between the stripes 
and superconductivity. The pairs give rise to a bond centered charge density periodicity 
of 4a along a wall which should give a Bragg peak at (j-, f-)- Other evidence of the pair 
structure of the domain wall should be seen in NMR measurements. 

The structure we have discussed differs from the one-electron mean field domain walls 
found in Hartree-Fock and Gutzwiller variational calculations in both its filling 
and its relationship to pairing. It also differs from the one-dimensional large U/t Hubbard 



model of a domain wall recently discussed by Nayak and Wilczek |I7 |. At J/t = 0.35, 
our system is well away from the phase separation regime, and we have no long-range 
intraplanar Coulomb interaction as in the frustrated phase separation domain wall model 



of Kivelson and Emery [|Tl|] . The singlet stripe phases discussed by Tsunetsugu [0,0 and 
the competing CDW-pairing phases in the 2- and 4-leg ladder systems recently discussed 
PJT0|] have some of the features found in this present study. In particular, they have the 
interplay of d x 2_ y 2 superconducting pairing and charge density wave order. Clearly, it will 
be important to understand if this is the behavior that underlies the structure of the domain 



walls in Lai. 6 _ a; Ndo.4Sr :r Cu04. 
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FIGURES 

FIG. 1. (a) Spin and hole structure suggested in Ref. [1] to account for neutron scattering 
experiments, (b) Hole density and spin moments for an 8 x 8 t-J model. The diameter of the gray 
holes is proportional to the hole density 1 — (n*), and the length of the arrows is proportional to 
(Sf), according to the scales shown, (c) For the same system, the exchange field strength — (SfSj) 
is given by the width of the lines according to the scale shown, when the holes (black dots) occupy 
their most probable configuration. In addition to showing all nearest neighbor exchange bonds, we 
show next nearest neighbor correlations about each hole if those correlations are antiferromagnetic, 
(Si ■ Sj) < 0. 

FIG. 2. (a) Hole density and spin moments for an 8 x 8 t-J model with a static potential 
AV = 0.05t applied to the sites in the shaded rectangles, (b) Charge structure factor S c (q) 
expected for a 2D system, obtained by periodically repeating the pattern shown in (a), (c) Spin 
structure factor S^q). Here we have measured I from the center of a domain wall. In both (b) 
and (c), the area of the circle is proportional to the structure factor, open circles denote negative 
values, and the axes are labeled in units of 1/a. 

FIG. 3. A single domain wall, modeled by a 4 x 18 system with left and right edge 7r-phase 
shifted antiferromagnetic fields h = 0.05i, with a static potential AV = 0.05i, and with an applied 
proximity affect field d = 0.03t. (a) Hole density and spin moments, (b) Pair field strengths on 
each nearest neighbor link, (c) The average pairing strength on the y links for the central 4 rungs, 
as a function of AV. 
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